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We study the superdeformed (SD) states and corresponding SD hypernuclei of Ar isotopes with the 
multidimensionally-constrained relativistic mean field (MDC-RMF) models which can accommodate 
various shape degree of freedom. We found that the density profiles of SD states in Ar isotopes show 
a strong localization with a ring structure near the surface, while the central part of the density is 
dilute showing a hole structure. This localization of SD density induces an appreciable deformation 
in the hyperon wave function and results in a large overlap between the core and the hyperon in 
the SD hypernuclei of Ar isotopes. Then the A separation energy of SD state becomes larger than 
that of normally deformed or spherical ground state. This feature is different from that found in 
other nuclei such as °S, °°Ni, and °°Zn in which the A separation energy of larger deformed state 
is smaller. In this context, the measurement of the A separation energy may provide an important 
information on the localization of the density profile of SD states. 


PACS numbers: 21.80.+a, 21.10.Dr, 21.60.Jz 


I. INTRODUCTION 


A A particle is free from the nucleon’s Pauli exclusion 
principle in nuclei due to the strangeness degree of free- 
dom. When a A particle is added to a nucleus, it lays 
deeply in the nucleus and many interesting phenomena 
were pointed out as the outcome. For instance, by an ad- 
dition of a A particle into the p-shell nuclei which have an 
a clustering structure, there is a dynamical contraction 
of nuclear size and stabilization in the binding energy [1— 
4|, which is called the glue like role of the A particle. 
However, the shrinkage effect is dependent on the states 
in the p-shell or sd-shell nuclei which have both a clus- 
tering states and shell-like compacting states. Namely, 
by the addition of a A particle, we have 20 % to 30 % 
nuclear shrinkage in the clustering states, while almost 
no shrinkage effect in shell-like states [5]. 

As an impurity in normal nuclei, it was pointed out 
that a hyperon also induces the change of the nuclear 
shape. The nuclear response to the addition of a A parti- 
cle can be studied by examining the A separation energy 
(Sa). The dependence of Sa on the nuclear shape is 
an interesting topic. For example, it has been pointed 
out that the A separation energies in clustering states of 
some nuclei, e.g., \?Be, 4°C, and 3'Ne, are smaller by 1 
to 3 MeV than those in shell-like compact states of cor- 
responding hypernuclei from few-body models [6, 7] and 
antisymmetrized molecular dynamics (AMD) model [8]. 


Recently, the self-consistent mean field models, both 
Skyrme Hartree-Fock (SHF) and relativistic mean field 
(RMF) models, have been applied to study the defor- 
mation of p- and sd-shell A hypernuclei [9-13]. In most 
of the cases, the deformation of core nuclei and the cor- 


responding hypernuclei are similar with the same sign. 
However, relativistic mean field models predict drastic 
changes of the deformation by the injection of a hyperon, 
e.g., in {8C and 49Si [10, 13]. Thus in general the polar- 
ization effect of hyperons is larger in RMF than Skyrme 
HF (SHF) models [11]. 


Beyond sd-shell nuclei, there are many interesting 
structure issues to study. One of them is that in some 
nuclei there exist superdeformed (SD) states which are 
characterized by the ratio 2:1 between the long and the 
short deformation axes in the coordinate space. In nuclei 
with A ~ 40, SD states have been observed experimen- 
tally [14-16] and the structure of these states has been 
extensively studied by mean field models [17-22]. 


A question related to the SD nuclei is, can the SD 
minima in normal nuclei persist after the injection of a 
hyperon? If the answer is yes, we have a further ques- 
tion, what happens about the nuclear responses in the 
normally deformed (ND) and SD states by the addition 
of a hyperon, say, a A? It is likely that the A separation 
energies in SD states are smaller than those in ND states 
because of smaller overlap between the A particle and SD 
core, as predicted for clustering states by the few-body or 
AMD models [6-8]. Especially, it has been pointed out 
that the A separation energies in SD states of Ca and 
Sc isotopes by the addition of a A particle are smaller 
than those of ND states by Isaka et al. [23]. However, 
it should be noticed that within the framework of mean 
field models, the wave functions of nuclear states might 
be localized [24, 25]. It would be interesting to investi- 
gate the energy gain by the addition of a A particle to 
such localized states by using mean field models. 


Among many mean field models, in this work we 
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adopt the multidimensionally-constrained relativistic 
mean field (MDC-RMF) models which can accommo- 
date various shape degrees of freedom [26] and have been 
used to study the shape of light hypernuclei [13] and the 
structure of heavy nuclei [27-29]. This paper is orga- 
nized as follows. In Sec. II, we describe briefly the MDC- 
RMF model and the properties of the A-N interaction 
employed here. In Sec. III, the results and discussions 
are presented. A summary is given in Sec. IV. 


II. THEORETICAL FRAMEWORK 


Recently multidimensionally-constrained relativistic 
mean field (MDC-RMF) models were developed in or- 
der to treat the various shape degrees of freedom in 
atomic nuclei. In the MDC-RMF models, the RMF 
equations are solved in an axially deformed harmonic 
oscillator (ADHO) basis [30]. The MDC-RMF models 
have been extended to including the hyperons [13], where 
the Dirac equations for the nucleons and hyperons and 
the Klein-Gordon equations for mesons are solved in the 
ADHO representation, while the Coulomb field is solved 
by the Green’s function method. Next we mention briefly 
the formalism of the MDC-RMF models, mostly those 
related to the A hyperon; the readers are referred to 
Ref. [26] for more details. 

In the MDC-RMF models, the RMF functional can 
be one of the following four forms: the meson exchange 
or point-coupling nucleon interactions combined with the 
nonlinear or density-dependent couplings. The deforma- 
tions 8), with even u can be considered simultaneously. 
In the present work we use the meson-exchange type in- 
teractions and only consider the axial and reflection sym- 
metric shapes. Thus the projection of the total angular 
momentum on the symmetric z-axis, Q, and the parity m 
are both conserved. 

The starting point of the meson-exchange RMF models 
for hypernuclei is the covariant Lagrangian density 


L= Llo + La, (1) 


where Lo is the standard RMF Lagrangian density for 
the nucleons and mesons [31-37] and La is that for the 
hyperon, 


Ly = Ya (itp — MA — Gord — Gury Wy) WA 


[DON 
—— ip 2 
F Ama Wao H Wa, ( ) 


where ma, is the mass of the A hyperon, goa and gwa are 
the coupling constants of the A hyperon with the scalar 
and vector meson fields, respectively. The last term rep- 
resents the tensor coupling between the A hyperon and 
the w field [38]. Quv is the field tensor of the w field 
defined as Quy = b pwy — O,w,. Couplings to the p me- 
son and the photon vanishes for A hyperons which are 
neutral and isoscalar. In the RMF model, the A — X° 


mixing, which was important in reproducing simultane- 
ously the binding energies of the s-shell hypernuclei (see, 
e.g., Refs. [39, 40]), is not included explicitly. 

Under the mean field approximation, the single particle 
Dirac equation for A hyperons reads, 


[a-p+B(ma+Sa)+Va t+ Tal dar = cipa, (8) 


with the scalar potential SA = gzao, the vector potential 
Va = gwaw, and the tensor potential, 


fona 
Tiers 
A IMA p (a 


p) w. (4) 


The mesons fields o, w, and p are obtained by solv- 
ing the Klein-Gordon equations with source terms. For 
example, the equation for the w meson reads, 


O Joh j (5) 


(-V?+m2)w = (gupp + Gwapnr) Oma PM 


where p% and pk are the isoscalar-vector densities for 
nucleons and hyperons and 


p = ið- (= vst) (6) 


i 


is the divergence of the current. 

In this work we use the parameter set PK1 [41] for 
the N-N interaction. The A-N effective interaction is set 
to be PK1-Y1 which was adjusted to reproduce the ob- 
served binding energies and spin-orbit splittings of the A 
hyper nuclei [42, 43]. In PK1-Y1 the coupling constants 
for the A-N interaction are determined by the relations 
goa = 0.580 go, gwa = 0.620 gu, and furan = —Gua. The 
spurious motion due to the breaking of the translational 
invariance is treated by including the microscopic center 
of mass correction Eem. = —(P?)/(2MA) in the total 
binding energy [41, 44, 45]. 

The pairing effects are included by the BCS approx- 
imation with a finite-range separable pairing interac- 
tion [46-48], 


Gi(R RPE =E, (7) 


where G is the pairing strength, R and F are the center 
of mass and relative coordinate between the paired par- 
ticles, respectively. P(7) is a Gaussian shaped function, 


1 =F? /a? 


AY) a aan” , 


(8) 
where a is the effective range of the pairing force. In 
Ref. [46] the strength and range of the separable force 
are adjusted to reproduce the momentum dependency of 
the pairing gap in the nuclear matter calculated from the 
Gogny forces. In this work we adopt the parameter set 
that mimics the Gogny force D1S: 


G = 728.0 MeV - fm?, a = 0.644 fm. (9) 
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FIG. 1: The energy of (a) ®Ar and Ñ Ar, (b) “Ar and 


3 Ar, and (c) *°Ar and 4'Ar as a function of the deformation 
parameter 32, respectively. The A hyperon is injected into the 
lowest s orbit and three p orbits, respectively. The Nilsson 
quantum numbers 2”[Nn3mz;] of the corresponding A orbits 
are marked. 


II. RESULTS AND DISCUSSIONS 


The calculated energy surfaces of several Ar isotopes 
(36 Ar, 33 Ar, and *°Ar) and the corresponding A hyper- 
nuclei are shown in Fig. 1. The A hyperon is injected into 


the lowest s orbit and three p orbits, respectively. The 
energy minimum of °° Ar is found at the oblate side with 
Bg ~ —0.2. The SD minimum is also found at 32 ~ 0.6. 
The ground state of 3° Ar is spherical as seen in Fig. 1, 
while the SD state locates at deformation 32 ~ 0.6 which 
is similar to that of 3°Ar. The energy surface of Ar is 
rather flat around {2 ~ 0 and the ground state is slightly 
oblate (see Table I). Around 82 ~ 0.5, there is a shallow 
minimum. 

The energy surfaces of hypernuclei rae (a hyperon 
is put in the 1sy/2 or Q*|Nngmj] = 1/2*[000] states) 
with A = 36, 38, and 40 essentially follow those of the 
core nuclei 4Ar. For instance, there are two minima for 
37 Ar, one at the oblate region and the other at the SD re- 
gion. The SD minima become rather shallow in the cases 
Ar and 4'Ar. Because of the small spin-orbit split- 
ting of hyperon-nucleon interaction, the 1p;/2 and 1p3/2 
states are almost degenerate in the limit of G2 = 0.0. The 
1/27 [110] and 3/27[101] states are splitted in energy by 
the 82 deformation in the similar way to the nuclear Nils- 
son levels while the 1/27 [101] state is almost degenerate 
with 3/27[101] state because of the very small spin-orbit 
splitting in the A potential. Around {2 ~ 0.15, there is 
a shoulder in the energy surface of °°Ar. With one A 
added, this shoulder persists; when the A occupies the 
1/27 [110] state, a very shallow minimum even develops 
around 62 ~ 0.15. In the following discussions, we will 
focus on the results in which the A hyperon occupies the 
1/2* [000] state. 

The calculated deformations of ground states and SD 
states of several Ar isotopes are listed in Table I together 
with the root mean square (r.m.s.) radii, the binding 
energies, and the A separation energy defined by 


Sa = BG Ar) = BA Ar). (10) 


From Table I one finds some familiar features concerning 
the impurity effects of the A hyperon. First, the shrink- 
age effect appears in all these nuclei studied here as seen 
from the values of the r.m.s. radii. This effect is quite 
small for the ground states: The radius of a hypernu- 
cleus is smaller by less than 0.02 fm than that of the 
core nucleus. Even for the SD states, the change of the 
radius does not exceed 1%. Second, with A increasing, 
the radius of A hypernuclei increases, so does the r.m.s 
radius of the A distribution. Third, the shape of a A 
hypernucleus always follows its core nucleus and the ad- 
dition of A does not change the sign of 62. Finally, once 
the core is deformed, the density distribution of A is also 
distorted and the deformation parameter Bə of a hyper- 
nucleus is always smaller than that of the corresponding 
core nucleus. 

In Table I, we find that the A separation energies of 
Ar isotopes at the ground states have smaller values com- 
pared with those of the SD states, which is opposite to the 
case of lighter A hypernuclei such as }°Be [6, 7]. The 1/27 
state in Be is much more deformed than the ground state 
3/21. When a A particle is added to these states, the A 
separation energy in °Be(1/27) + A state is smaller than 
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TABLE I: The quadrupole deformation parameters, root mean square (r.m.s.) radii, binding energies, and A separation energies 
for the ground states and SD states (labelled with asterisks) in several Ar isotopes and the corresponding A hypernuclei. The 
experimental values of binding energies are taken from Ref. [49-51]. 


Nucleus Quadrupole deformation parameters R.m.s. radii (fm) Energies (MeV) 

b2 Bn Bp Ba Tm TA E Eexp Sa 
Ar 0.212 —0.208 —0.215 3.238 —303.802  —306.716 
37 Ar —0.204 —0.205 —0.211 —0.057 3.220 2.644 —321.979 18.177 
38 Ar 0.000 0.000 0.000 3.281 —326.455  —327.343 
3 Ar 0.000 0.000 0.000 0.000 3.265 2.672 —344.896 18.441 
40 Ar —0.123 -—0.112 —0.137 3.338 —342.613 —343.810 
“Ar —0.117 —0.109  -—0.132 —0.037 3.321 2.698 —361.398 18.785 
35 Are 0.620 0.610 0.630 3.346 —296.670 
3T Ar” 0.597 0.599 0.619 0.172 3.319 2.626 —315.194 18.524 
38 Ar* 0.602 0.597 0.609 3.403 —317.448 
9° Ar*® 0.589 0.596 0.603 0.187 3.378 2.659 —336.306 18.858 
40 Ar* 0.499 0.472 0.533 3.430 —336.852 
TAr" 0.491 0.474 0.530 0.161 3.409 2.685 —355.922 19.070 
TABLE II: Calculated deformation parameter 82, the over- delormaribw patdnmieter Ha, the binding energy Æ, and A 


lap Ioverlap defined in Eq. (11), the binding energy E, and A 
separation energy Sa of some A hypernuclei. For comparison, 
the binding energy of the core nucleus Ecore is also given. The 
energies are in MeV. 


Nucleus b2 Toveilap Erot Fore SA 

3T Ay —0.204 0.1352 —321.979 —303.802 18.177 
37 Ar* 0.597 0.1370 -—315.194 -—296.670 18.524 
39 Ar 0.000 0.1360 -—344.896 -—326.455 18.441 
39 Ar* 0.589 0.1378 —336.306 —317.448 18.858 
“Ar —0.117 0.1357 —361.398 -—342.613 18.785 
aT Ar* 0.491 0.1378 —355.922 -—336.852 19.070 
“Ca 0.00 0.1361 361.422 —342.869 18.553 
ICat 0.70 0.1393 —350.559 —331.317 19.242 
pas) 0.26 0.1376 —285.095 —267.002 18.093 
33 9* 0.97 0.1243 —274.315 —257.951 16.364 
Ni 0.00 0.1461 —506.665 —484.759 21.906 
STNi* 0.40 0.1415  —498.610 —477.892 20.718 
“Zn 0.22 0.1438 —534.565 —512.924 21.641 
ĉl Zn* 0.62 0.1415 —527.168 —506.238 20.930 


that in °Be(3/2;) +A state. This has been explained as 
that the overlap between nucleons and the A hyperon in 
the 1/2} state is smaller than that in the ground state. 
Next we will show that in Ar isotopes, however, the over- 
lap between the nucleons and the A hyperon in the SD 
state is larger than that in the ground state, which results 
in a larger A separation energy in the SD state. 

To clarify the physical mechanism of the larger A sep- 
aration energy in the SD state, we calculate the overlap 
Ioverlap Of the core and the hyperon, 


Ioverlap = [Pool z)pa(r,z)rdrdz, r= v/a? + y?. 
(11) 


where the core and the hyperon densities are denoted by 
Pcore(r, z) and pa(r, z), respectively. 
In Table II, the overlaps overlap are listed together the 


separation energy S, of some A hypernuclei. For ‘ae? 
one finds that the overlap Joverlap in the SD state is always 
larger than that in the ground state. This explains well 
that the A separation energy in the SD state is larger 
than that of the ground state. In the following we will 
examine the density distributions of the core nucleus and 
the A hyperon and show that this new feature actually 
stems from a localized density of nucleons in the core 
nucleus and a stretched density of the A hyperon in the 
SD state. 


We illustrate the calculated two-dimensional density 
distributions in the r-z plane (note that the z-axis is the 
symmetric one) for the ground and SD states in 36:38-40 Ar 
and 4"9-4! Ar in Fig. 2. The ground states of 3°Ar and 
its hyper-counterpart 3”Ar are clearly oblate. On the 
other hand, it is remarkable that the SD density dis- 
tribution shows a clear localization feature with a ring 
shape near the surface and a hole structure with lower 
density at the center. Note that similar localization fea- 
ture has been found and studied in details for the ground 
state of 2°Ne [24]. Furthermore, we see that as increasing 
number of neutrons, this localization with a ring shape 
of density profile becomes less pronounced as shown in 
Fig. 2. This behavior can also be deduced in {2 values in 
Table I where one finds that, as the number of neutrons 
increases, the deformation 32 becomes smaller. When a 
A particle is added to 36:38-40 Ar, the density profiles are 
almost similar with those of corresponding normal nu- 
clei and the localization effect is still there. And also, as 
increasing the number of neutrons, the weakness of local- 
ization with a ring shape of density profile in hypernuclei 
is similar with those of the corresponding normal nuclei. 


In Fig. 3 we show the two-dimensional density distri- 
butions in the r-z plane for the A hyperon in the ground 
and SD states of %’Ar. It is seen that in the ground of 
4’ Ar, the A hyperon distributes almost isotropically; this 
is consistent with the small 6, value given in Table I. 
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FIG. 2: Two-dimensional density distributions in the r-z 
plane (r = yx? + y? and z-axis is the symmetric one) for (a) 
36 Ar and 4" Ar and (b) *8Ar and $° Ar and (c) Ar and 4! Ar 
The asterisks are marked for SD states. 
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FIG. 3: Two-dimensional density distributions in the r-z 
plane (r = \/x? + y? and z-axis is the symmetric one) for the 
A hyperon in the ground state 4’ Ar and the SD state 3’ Ar*. 


TABLE III: The probabilities of A particle wave function in 


different Nilsson orbits for 4’ Ar in ground state and SD state, 

respectively. 

Nilsson levels GS SD 
1/27 [000 0.955 0.950 
1/27 [220 0.006 0.038 
1/27 [200 0.034 0.006 
1/27 [440 0.000 0.002 


Therefore in the ground state of 3 Ar, the overlap be- 
tween the core and the hyperon wave function is small 
because the core has an oblate shape while the hyperon 
is almost spherical. However, in the SD state of STAT, 
the density distribution of the A hyperon is stretched 
substantially along the z-axis and the quadrupole defor- 
mation parameter 6, = 0.172 (see Table I). The proba- 
bilities of the hyperon wave function in different Nilsson 
orbits are given in Table III for the ground and SD states 
in 37 Ar. It is seen that the hyperon wave function embed- 
ded in the SD state has larger high-N components (even 
4ħw excitation) than in the ground state. These higher- 
N components with m; = 0 corresponds to larger prolate 
shape [52, 53] and results in a larger overlap between the 
hyperon and the core in 37 Ar*. Since the SD density dis- 
tribution shows a localization feature with a ring shape 
near the surface and a hole structure with low density 
at the center and the radial coordinate r is included in 
the integrand, i.e., TPcore(r, Z)pa(r, z) in Eq. (11), it is 
natural that the overlap [overlap is larger in the SD state 
than that in the ground state of 3’ Ar. 

To clarify the relation between the large A-separation 
energy in the SD state and the localization feature of 
the core nucleus, we investigate nuclei which have SD 
minima in their potential energy surfaces such as °7S, 
49Ca, 56Ni, and ©°Zn and the corresponding A hypernu- 
clei. Among these nuclei, as shown in Fig. 4, for Ca 
the SD state has a localization in its density profile with 
a ring shape and in other SD states there are no local- 
ization effects. Then it is expected that the localization 
in the SD state causes again a large overlap between the 
core and the hyperon wave function and, as a result, the 
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FIG. 4: Two-dimensional density distributions in the r-z 
plane (r = \/x? + y? and z-axis is the symmetric one) in the 
ground state and the SD state of *°Ca. 


A separation energy of the SD state becomes larger than 
that of the ground state. On the other hand, it is likely 
that since other isotopes, 32S, °°Ni, and ©°Zn have non- 
localized core densities, the A separation energies in SD 
states should be smaller than those in the ground states. 
In Table II, we give the overlaps Joverlap together the de- 
formation parameter (2, the binding energy E, and A 
separation energy Sa of these A hypernuclei. For a Ca, 
the overlap Ioverlap in the SD state is indeed larger than 
that in the ground state, so is Sa. However, for 33S, 4’ Ni, 
and 8'Zn, the overlap between the hyperon and the SD 
core is smaller than that of the ground state which is 
either normally deformed or spherical. Then the A sep- 
aration energy is larger for the smaller deformed state 
than for the SD state. 

From the above discussions, we conclude that the local- 
ization in the SD density tends to derive an appreciable 
deformation in the hyperon wave function (see Fig. 3) 
and a larger overlap between the core and the hyperon 
which in turn results in a larger A separation energy. 
However, if in the SD state there is no localization in the 
core density, the overlap between the core and the hy- 
peron is smaller in the SD state than that in the ground 
state, which results in a larger A separation energy in 
the latter. These conclusions have been partly demon- 
strated in recent AMD model calculations by Isaka et 
al.: The density for the SD state in 3ĉAr does not show 


1] T. Motoba, H. Bando, and K. Ikeda, Prog. Theor. Phys. 
70, 189 (1983). 

2| T. Motoba, H. Bando, K. Ikeda, and T. Yamada, Prog. 
Theo. Phys. Supp. 81, 42 (1985). 

3] E. Hiyama, M. Kamimura, T. Motoba, T. Yamada, and 
Y. Yamamoto, Phys. Rev. C 53, 2075 (1996). 

4) E. Hiyama, M. Kamimura, K. Miyazaki, and T. Motoba, 
Phys. Rev. C 59, 2351 (1999). 

5| E. Hiyama, M. Kamimura, T. Motoba, T. Yamada, and 
Y. Yamamoto, Prog. Theo. Phys. 97, 881 (1997). 

6] E. Hiyama and Y. Yamamoto, Prog. Theo. Phys. 128, 
105 (2012). 


7| Y. Zhang, E. Hiyama, and Y. Yamamoto, Nucl. Phys. 


the localization feature and the A separation energy of 
the ground state is larger than that of the SD state in 
4’ Ar [54] which is contradictory to the present results 
from RMF models. It would be very interesting to study 
the localization feature of different mean field models for 
nuclear SD states. 


IV. SUMMARY 


We studied the SD states and corresponding SD hyper- 
nuclei of Ar isotopes by using the RMF model. We found 
that the density profiles of SD states in Ar isotopes show 
a strong localization with a ring structure near the sur- 
face, while the central part of the density is dilute show- 
ing a hole structure. The strong localization of the SD 
state is also found in 4°Ca. This localization feature of 
the SD density induces an appreciable deformation in the 
hyperon wave function of the SD hypernuclei. Then the 
A separation energy of SD state becomes larger than that 
of normally deformed or spherical ground state. This fea- 
ture is different from that found in other nuclei such as 
$38, 9’Ni and $'Zn in which the A separation energy of 
the SD state is smaller. In this context, the measurement 
of the A separation energy would provide an important 
information on the localization of the density profile of 
SD states. 
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